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ABSTRACT 

We measure the average gravitational shear profile of 6 massive clusters (M v ; r ~ 10 15 M Q ) at z = 0.3 out to a 
radius ~ 2h Mpc. The measurements are fitted to a generalized NFW-like halo model p(r) with an arbitrary 
r — > slope a. The data are well fitted by such a model with a central cusp with a ~ 0.9- 1.6 (68% confidence 
interval). For the standard-NFW case a = 1.0, we find a concentration parameter c v ; r that is consistent with 
recent predictions from high-resolution CDM N-body simulations. Our data are also well fitted by an isothermal 
sphere model with a softened core. For this model, our ler upper limit for the core radius corresponds to a limit 
cr* < 0. lcm 2 g _1 on the elastic collision cross-section in a self-interacting dark matter model. 

Subject headings: dark matter — gravitational lensing — Galaxies: clusters 
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1. INTRODUCTION 

Although the cold dark matter (CDM) model of structure 
formation has been successful in explaining many observable 
properties of the universe, there remain some notable discrep- 
ancies with observations. In particular, numerical simulations 
of dark matter haloes in CDM universes appear to overpro- 
duce the abundance of dark matter sub-haloes corresponding 
to dwarf galaxies by at least an order of magnitude (Klypin et 
al. 1999; Moore et al. 1999). Also, simulations indicate a trans- 
fer of angular momentum from the baryonic matter component 
to the non-baryonic dark matter during the assembly of disk 
galaxies, which leads to a conflict with the observed angular 
momentum properties of disk galaxies (e.g., Navarro & Stein- 
metz 1997,2000; Sommer-Larsen, Gotz, & Portinari 2002). 

Finally, the simulations predict that dark matter halos in a 
wide mass range from dwarf galaxies to massive clusters fol- 
low an universal density profile with a central cusp p(r) oc r~ a 
with a in the range 1.0 < a < 1.5 (Navarro, Frenk, & White 
1996, 1997; hereafter collectively NFW; Moore et al. 1999). 
This appears to be contradicted by some observations that in- 
dicate an almost constant-density core rather than a cusp, both 
on galaxy scales (e.g., Dalcanton & Bernstein 2000; Salucci & 
Burkert 2000; de Blok et al. 2001) and on cluster scales (Tyson, 
Kochanski, & Dell'Antonio 1998), although the interpretation 
of these results is controversial (e.g., Broadhurst et al. 2000; 
Shapiro & Iliev 2000; Czoske et al. 2002). This has led to sug- 
gestions that the dark matter properties may deviate from stan- 
dard CDM, e.g., the dark matter could be warm (Colin, Avila- 
Reese, & Valenzuela 2000; Sommer-Larsen & Dolgov 2001), 
repulsive (Goodman 2000), fluid (Peebles 2000), fuzzy (Hu, 
Barkana, & Gruzinov 2000), decaying (Cen 2001), annihilating 
(Kaplinghat, Knox, & Turner 2000), self-interacting (Spergel 
& Steinhardt 2000; Yoshida et al. 2000, Dave et al. 2001) or 
both warm and self-interacting (Hannestad & Scherrer 2000). 
Alternatively, it has been suggested that stellar feedback from 
the first generation of stars formed in galaxies was so efficient 
that the remaining gas was expelled on a timescale comparable 
to, or less than, the local dynamical timescale. The dark matter 
subsequently adjusted to form an approximately constant den- 
sity core (e.g., Gelato & Sommer-Larsen 1999). This is how- 



ever unlikely to affect cluster cusps. 

It is dangerous to draw wide-ranging conclusions about halo 
density profiles from a single object, even when based on high- 
quality data from strong gravitational lensing. Clearly, a larger 
number of clusters should be investigated using a range of 
methods. Other studies have sought to put constraints on a 
from the observed abundances of strongly lensed arcs produced 
by clusters (Meneghetti et al. 2001; Molikawa & Hattori 2001; 
Oguri, Taruya, & Suto 2001), and find best fit values of a simi- 
lar to CDM predictions. 

Previous weak lensing measurements of cluster mass profiles 
have usually only considered two models, a singular isothermal 
sphere (SIS) model with p oc r" 2 and the NFW model. Typ- 
ically, it is found that the NFW model provides a marginally 
better fit to the data, but usually not at a very high significance 
(Clowe et al. 2000; Clowe & Schneider 2001; Hoekstra et al. 
2002). Many of the clusters studied in this way show evidence 
for substructure, indicating that they may not be in dynamical 
equilibrium. 

Here, we consider a density profile which is a generalized 
version of the NFW model, as described in §2. The data 
we present here, detailed in §3, consists of measurements of 
gravitational shear over a wide range of radii around 6 clus- 
ters at z — 0.30. In §4 we present the results of our model 
fits, and in §5 we discuss the implications for cluster-scale 
dark matter halos and dark matter physics. We consider both 
a spatially-flat ACDM model with (fi ,A ) = (0.3,0.7) and a 
SCDM (Einstein-de Sitter) cosmology with (Sl , A ) = (1,0), 
and we use Ho = 100/zkms _1 Mpc . 

2. HALO MODEL 

Here we consider a generalized halo density profile on the 
form (Zhao 1996; Jing & Suto 2000) 



ScPc 



(1) 



(r/r s ) Q (l + (r/r. s )) 3 - Q ' 

which corresponds to the standard-NFW model for a = 1 , but 
has an arbitrary inner slope. Moore et al. (1999) favor an inner 
slope a = 1 .5 and adopt a model with a slightly different func- 
tional form than we consider here, but with the same asymptotic 
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behavior as an a = 1 .5 model in equation Q for large and small 
scales. 

NFW originally defined a concentration parameter c by c = 
7"20oA.s, where /5( r 20o) = 200p c and p c is the critical density. For 
the model above, the characteristic density S c is 



5c = 



200 



-i -l 



x 2 (cxT a (l+cx) a - 3 dx 



(2) 



Keeton & Madau (2001) define a slightly different concentra- 
tion parameter c_2 = r v i r /r_2, where r v i r is the virial radius of the 
halo and r_2 = (2-a)r s is the radius at which din p/ din r = -2. 
At the redshift of the clusters we examine (z = 0.3), N-body sim- 
ulations of clusters in a ACDM universe indicate that the virial 
radius is related to r 2 oo by r vil - = 1 . 194r 2 oo+0.003Mpc (M. Gotz, 
2002, private communication). For the SCDM model we make 
the approximation r v h- ~ r2oo, since p(r v i r ) = 178p f ~ 200p c . 

Bullock et al. (2001) use high-resolution N-body simula- 
tions to predict the mass- and redshift-dependence of the con- 
centration parameter c vlr = r w - u /r s for NFW-type halos. Sim- 
ulated dark matter halos show a significant dispersion about 
the median value of the concentration parameter at a given 
mass (Jing 2000), and Bullock et al. (2001) predict a scatter 
of A(logc vil ) = 0.18. Bullock et al. (2001) found the redshift- 
dependence of this median value to be c v ; r oc (1 + z) _1 . For 
very massive cluster halos (M v ; r ~ 10 15 M Q ), the concentration 
is low, and r„ is of order several hundred kpc, meaning that the 
difference between various values of a is evident even at cluster 
radii associated with the weak lensing regime. 

3. LENSING DATA 

The cluster sample used here was selected from a larger sam- 
ple of 38 clusters targeted for weak lensing mass measurements. 
The data set, and the methods for data reduction and shear es- 
timation have been described elsewhere (Dahle et al. 2002). In 
this letter we focus on the density profiles of a subset of six clus- 
ters of galaxies at z = 0.3. These are the most distant clusters in 
the sample of Dahle et al. (2002) for which wide-field imaging 
data exist, allowing us to probe the cluster mass profile beyond 
the virial radius. 

The data were obtained with the UH8K mosaic CCD cam- 
era at the 2.24 m University of Hawaii telescope at Mauna Kea 
Observatory (see Dahle et al. 2002 for details). The clusters 
were selected from the X-ray luminous cluster samples of Briel 
& Henry (1993) and Ebeling et al. (1996), and all have X-ray 
luminosities L x ,o.5-2.5keV > 4 x 10 44 erg s" 1 . The virial masses 
M v h of half of the clusters have been estimated by Irgens et al. 
(2002), who find values in the range 1.1 x 10 15 < hM vil /M Q < 
3.3 x 10 15 . 

The gravitational shear 7 and convergence k produced by the 
halo model given in equation Q cannot be expressed in terms 
of simple functions, but we have derived approximate relations 
that are accurate to about 2% for the part of parameter space that 
we are able to probe with our data. The observable tangential 
distortion gj is related to the tangential shear and convergence 
by gr = 7r/(l ~~ K ) at all radii larger than the Einstein radius r^. 

As first suggested by Kaiser & Squires (1993), the signal- 
to-noise ratio of weak lensing observations can be boosted by 
"stacking" the weak shear measurements of a number of sim- 
ilar objects. Here, we stack our tangential distortion measure- 
ments for the clusters A959, A1351, A1705, A1722, A1995 and 
A2537, which all lie in the redshift interval 0.285 < z < 0.325. 



For each cluster, the average tangential distortion (gr) vs. ra- 
dius was measured in 24 radial bins, assuming the peak in the 
reconstructed projected mass density (see Dahle et al. 2002) to 
be the cluster center. 
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FIG. 1 . — The points with error bars show the averaged gravitational lens- 
ing distortion of background galaxies as a function of radius. The solid line 
indicates the best fit generalized-NFW model and the dashed line indicates the 
best fit NSIS model. Corrections for cluster member contamination have been 
applied to the model curves. The physical radii plotted here are calculated for 
the SCDM cosmology and would be 15% larger for the ACDM cosmology. 

We re-scale our clusters to the same mass (represented by 
a halo with velocity dispersion a v = 1 150km s -1 , which is the 
average value for our clusters) using a SIS model fit to the 
observed tangential distortion. In doing this, we neglect the 
mass-dependence of the concentration parameter and the c v ; r oc 
(z+ 1) -1 redshift dependence predicted by Bullock et al. (2001). 
These are both only minor corrections, since the clusters we 
consider here have a very small spread in redshift and modest 
scatter in mass, and at this point we do not want to impose addi- 
tional constraints on our model parameters by assuming scaling 
laws like those predicted by NFW or by Bullock et al. (2001). 
We also scale the clusters to a uniform lens redshift z = 0.3 by 
taking into account the variation in lensing strength and angular 
diameter with redshift. 

The faint galaxy catalogs which were used to estimate (gr) 
were generated by selecting galaxies that were detected at low 
significance v (6 < v < 100). These galaxies have typical red- 
shifts Zb g ~ 0.9 and typical magnitudes V = 24, / = 23. At small 
projected radii from the cluster center there will be a signifi- 
cant fraction of cluster galaxies in the faint galaxy catalogs, and 
the shear estimates will be systematically biased towards lower 
values due to this contamination. The importance of this effect 
can be estimated by assuming that there is negligible contribu- 
tion by cluster galaxies at the edges of the UH8K fields (i.e., 
beyond ~ 1.5/T 1 Mpc). Thus, the true density of background 
galaxies can be estimated, and the cluster contribution to the 
faint galaxy density can be calculated in azimuthally averaged 
bins as a function of radius from the cluster center. By cal- 
culating the average contamination for the six clusters that go 
into the stack, we can significantly reduce the effects of cluster 
richness variations, substructure and background galaxy den- 
sity fluctuations. 

4. RESULTS 

The average tangential distortion profile in a SCDM universe 
is shown in Figure ^ along with the best fit generalized-NFW 
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model curve. The 68% confidence interval for each parameter 
is given in Table ^for the two cosmologies we consider. The 
observed profile is well represented by equation Q for both the 
SCDM and ACDM case, with best fit x 2 /d.o.f. = 23.66/21 and 
22.78/21, respectively. 
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FIG. 2. — The contours indicate the 68% and 95% joint confidence intervals 
on c v ; r and a, derived from variations in \" w i'h respect to the minimum value. 
Also shown is the mean value and scatter in the concentration parameter for a 
M vir = 2 X 10 15 M Q NFW halo predicted by Bullock et al. (2001) for a ACDM 
cosmology. The dashed lines indicate lines along which the two parameters 
are degenerate. 

A massive central galaxy can significantly contribute to the 
projected mass density at small radii. Williams, Navarro & 
Bartelmann (1999) found that strong lensing constraints on 
clusters with a v ~ 1500- 2000km s~ are consistent with a 
standard-NFW profile, but a v ~ 1000km s _I clusters are incon- 
sistent with this model, unless the central surface density is sig- 
nificantly enhanced by the central galaxy. We check the robust- 
ness of our results to the inclusion of such an additional mass 
component by adopting a combined model which is a superpo- 
sition of the dark matter profile in equation Q and a SIS mass 
profile with a v = 300km s . This velocity dispersion is a typi- 
cal value measured for brightest cluster galaxies (e.g., Oegerle 
& Hoessel 1991) and is consistent with the extra central mass 
of ~ 3 x lO I2 /2 -1 M required to produce the arcs studied by 
Williams etal. (1999). 
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FIG. 3. — Constraints on the core radius r c for a NSIS profile (solid line: 
SCDM; dashed line: ACDM). 

As shown in Tabled the inclusion of an additional SIS mass 
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profile did not significantly alter the best fit parameter values, 
and the quality of the fits was also not strongly affected. Not 
surprisingly, it is seen that with the inclusion of an extra mass 
component with a steep r~ 2 slope in the center, the best-fit inner 
slope of our generalized NFW model needs to become slightly 
shallower to match the observed shear. Figure [2] shows 2- 
parameter joint confidence regions for c wu and a for the ACDM 
cosmology. For the a = 1 case, the confidence intervals for c vrr 
are found to be consistent with the range of concentration pa- 
rameters for M v ir = 2 x 10 15 )i~ 1 Mq halos predicted from N-body 
simulations by Bullock et al. (2001). 

Table 1 
Best fit model parameters. 



SCDM(n m = l, 
Generalized-NFW 
Parameter Best fit value 68% conf. int. 
c vir 1.3 0.2-4.2 

a 1.5 0.9-1.6 

r vir 1.34 1.28-1.49 

ACDM (U m = 0.3, Q A 
Generalized-NFW 
Parameter Best fit value 68% conf. int. 
c vir 1.7 0.4-2.6 

a 1.4 1.3-1.6 

r vir 1.97 1.91-2.03 



n A = o) 

Generalized-NFW + SIS 
Best fit value 68% conf. int. 

1.3 0.3-2.8 

1.4 1.1-1.6 
1.34 1.28-1.43 

= 0.7) 

Generalized-NFW + SIS 
Best fit value 68% conf. int. 
2.1 0.2-4.8 
1.3 0.9-1.6 
1.91 1.79-2.10 



We also fitted our data to a non-singular isothermal sphere 
(NSIS) model on the form p(r) = a 2 /[2tt G(r 2 + r 2 )], where r c 
is the core radius. Our data were also well described by this 
model, with a best fit x 2 /d.o.f. = 23.5/22 for both SCDM 
and ACDM. Our best fit values were r c = 58/T 1 kpc and r c = 
66/T 1 kpc for SCDM and ACDM, respectively. The best fit 
model curve is plotted in Figure [2 and the resulting 68% and 
95% confidence intervals for r c are shown in Figure[3] 

5. DISCUSSION 

We find that the average density profile of 6 massive clus- 
ters is well fitted by a generalized-NFW model with central 
slope a ~ 0.9- 1.6. This (steep) slope fits well with predic- 
tions from collisionless CDM simulations. The best fit slope of 
a = 1.3— 1.5 is slightly steeper than the canonical NFW value, 
although a = 1 is only disfavored at less than 90% confidence 
for any of the models. For ACDM, a < 0.5 can be excluded 
with 95% confidence. The inclusion of an additional SIS mass 
profile corresponding to a massive central galaxy only resulted 
in minor changes for the best fit parameters. This indicates 
that any additional mass component associated with the central 
galaxy does not make a significant contribution to the total pro- 
jected density on the scales that we probe in our weak lensing 
study. 

Our constraints on a are compatible with recent strong lens- 
ing studies of massive clusters: Smith et al. (2001) estimate 
a = 1.3 for A383 and Gavazzi et al. (2002) find 0.7 < a < 1.2 
for MS2 137-23 (see, however, Sand, Treu & Ellis 2002, who 
find a best fit a = 0.35 and a < 0.9 at 99% CL for the same 
cluster). Our results are also consistent with recent X-ray clus- 
ter studies based on Chandra data: Arabadjis, Bautz, & Garmire 

(2002) find that the mass profile of CL 1358+6245 is well fit 
by a standard-NFW a = 1 profile, and Lewis, Buote, & Stocke 

(2003) find a = 1.19 ± 0.04 for the central slope of A2029. 
We find that an isothermal sphere model with a softened core 

can also produce a good fit to our data. This indicates that our 
weak lensing data is unable to discriminate between an outer 
slope of r~ 2 and r" 3 . A pure SIS model is however strongly 
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excluded at the 99.9% CL (see Fig. 0. Yoshida et al. (2000) 
use simulations to study the effect of self-interacting dark mat- 
ter on the structure of cluster halos. They find that even a 
small elastic collision cross-section of cr* = 0.1cm 2 g _1 will sig- 
nificantly affect the central density profile. Meneghetti et al. 
(2001) study the strong lensing properties of the cluster mod- 
els of Yoshida et al. (2000), and find that even a cross-section 
as small as 0.1cm 2 g _1 is incompatible with the observed abun- 
dances of strongly lensed radial and tangential arcs in clusters. 
Their study is based on a simulated cluster at z = 0.278 of final 
virial mass 7.4 x 10 14 /r'Mo in a ACDM universe. This cluster 
is at a similar redshift as the clusters in our data set, but our clus- 
ters are on average about twice as massive. From the r c values 
tabulated by Meneghetti et al. (2001), we find that the predicted 
core radius for cr* = 0. lcm 2 g _1 of r c = 80/T 1 kpc is similar to our 
lcr (and 2er) upper limit for the core radius shown in Figure [3] 
Here, we have taken into account a factor y3 difference in the 
definition of r c and a M 1 / 3 scaling of lengths. A cross-section 
of lcm 2 g _1 would produce a core radius of 180/T 1 kpc, which 
is strongly excluded by our data. In the small mean free path 
limit of "fluid" dark matter, a core-collapse producing a SIS- 



type profile is expected. This is strongly excluded by our data, 
and the fluid limit is also excluded by the observed ellipticities 
of clusters (Miralda-Escude 2002). 

Dave et al. (2001) show that cross-sections of order 5cm 2 g _1 
are needed to produce a good fit to the large apparent cores 
in dwarf galaxies. By introducing a velocity-dependent in- 
teraction cross-section for the dark matter particles, it may 
still be possible to reconcile our results with a self-interacting 
dark matter model that can produce constant-density cores in 
dwarf galaxies. Even if such a model may still be allowed, 
its velocity-dependence is already strongly confined by other 
astrophysical constraints, including the demographics of su- 
permassive black holes in galaxy nuclei (Hennawi & Ostriker 
2001). 

It should be noted that even though dark matter models which 
produce large constant-density cores on clusters scales appear 
to be ruled out, many of the non-standard dark matter models 
are easily compatible with cusped cluster halos. One example is 
warm dark matter where the thermal velocity of the dark matter 
particles can erase the cusps of dwarf galaxy halos, but almost 
not affect more massive halos at all. 
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